Response to tilted magnetic fields in Bi2Sr2CaCu20s with columnar defects : Evidence 

for transverse Meissner eff'ect. 
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The transverse Meissner effect (TME) in the highly layered superconductor Bi2Sr2CaCu20s+y 
with columnar defects is investigated by transport measurements. We present detailed evidence for 
the persistence of the Bose-glass phase when H is tilted at an angle 6 < Oc (T) away from the column 
direction: (i) the variable- range vortex hopping process for low currents crosses over to the half- 
loops regime for high currents; (ii) in both regimes near Oc{T) the energy barriers vanish linearly 
with tan6 ; (iii) the transition temperature is governed by Tbg{0) — TBciO) ~ |tan0|^'''^^ with 
v±_ = 1.0 ± 0.1. Furthermore, above the transition as 9 ^ Ot , moving kink chains consistent with a 
commensurate-incommensurate transition scenario are observed. These results thereby clearly show 
the existence of the TME for 6 < OriT) . 

PACS numbers; 74.60Ge , 74.72Hs 



The mixed state in high-Tp superconductors has been 
of considerable interest over the past ten years. The rea- 
son is that combined effects between temperature, elas- 
ticity, intervortex interactions, and disorder yield a large 
variety of new thermodynamic vortex phases ^ . In par- 
ticular, the theory of vortex pinning by correlated dis- 
order such as twin planes or artificial columnar defects 
has been considered by Nelson and Vinokur (NV) and 
Hwa et al In the case of parallel columns, NV have 
shown that if the applied magnetic field H is aligned with 
the columnar defects, the low temperature physics is sim- 
ilar to that of the Bose glass (BG) [Q, with the flux lines 
strongly localized in the tracks leading to zero resistivity. 
When H is tilted at an angle 6 away from the column 
direction, the BG phase with perfect alignment of the in- 
ternal flux density B parallel to the columns is predicted 
to be stable up to a critical transverse field Hj_c produc- 
ing the so-called transverse Meissner effect (TME) [^1). 
For 9 > 9c = arctan(/ioi?±c/S), the linear resistivity is 
finite resulting from the appearance of kink chains along 
the transverse direction as discussed in Ref . |^ . Finally, 
above a still larger angle, the kinked vortex structures 
disappear and B becomes collinear with H. Similar sce- 
narios have also been developed for pinning of vortices 
by the twin boundaries H] and by the layered structure 
of the coumpound itself |0| . 

Recently, the TME in untwinncd single crystals of 
Y Ba2CuzO'j (YBCO) with columnar defects, has been 
observed using an array of Hall sensors js). Previ- 
ous magnetization or magnetic-torque measurements in 
anisotropic three-dimensional superconductors, such as 
YBCO, gave support to the presence of vortex lock- 
in phenomena, due to pinning by the twin boundaries 
or by the interlayers between the Cu-0 planes [ p^ . 
In the case of highly layered superconductors, such as 
Bi2Sr2CaCu20s+y (BSCCO), although the lock-in tran- 
sition was observed for H tilted away from the layers [O , 



the existence of the TME due to the pinning by columnar 
defects remained an open question. 

In this letter, we present measurements of the electri- 
cal properties near the glass-liquid transition in BSCCO 
single crystals with parallel columnar defects, as a func- 
tion of T and 9 for filling factors / < 1. We find that 
when 9 is less than a critical angle 9c{T), the resistiv- 
ity p(J) vanishes for low currents. A detailed quantita- 
tive analysis of p( J) indicates that the creep proceeds via 
variable-range vortex hopping (VRH) at low currents due 
to some on-site disorder [|2[, crossing over to the half- 
loop regime at high currents. For 9 > 9c(T), the very 
signature of a kinked vortex structure, consistent with 
a commensurate-incommensurate (CIC) transition sce- 
nario in (1+1) dimensions Q, is deduced from the mea- 
surement of the critical behavior of the linear resistivity. 
All these results clearly demonstrate that when the field 
is tilted at 9 < 9c away from the defects, the flux lines 
remain localized on columnar defects, and hence, the BG 
exhibits a TME. 

The BSSCO single crystal was grown by a self-flux 
technique, as described elsewhere |l^. The crystal of 
1 X 1 X 0.030 mm^ size with the c-axis along the short- 
est dimension has a Tc of 89 X, a transition width of 
^ 1 K, and was irradiated along its c-axis with 5.8 GeV 
Pb ions to a dose corresponding to 5$ = 1.5 T at 
the Grand Accelerateur National dTons Lourds (Caen, 
France) . The crystal was mounted onto a rotatable sam- 
ple holder with an angular resolution better than 0.1° in 
a cryostat with a 9 T magnet. Isothermal I — V curves 
were recorded using a dc four-probe method with a sen- 
sitivity of ~ 10"^" V and a temperature stability better 
than 5 mK. H was aligned with the tracks using the 
well-known dip feature occurring in dissipation process 
for 9 = 0°, and was tilted away from the column direc- 
tion for / = ^qH\\/B^ fixed. We present in this paper 
the data obtained for / = 1/3. 
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Fig. (1) shows a typical log-log plot oiV/I — I curves 
obtained varying 9 ior T < Tbg{^) fixed. We observe 
a well-defined angular crossover at an angle dc- For low 
angles {9 < 9c), V/I~ I data suggest a glassy vortex- like 
state resulting from the persistence of the BG phase pre- 
viously observed in this crystal for 9 — 0° [|l2|. On the 
contrary, upper curves {9 > 9c) display an Ohmic regime, 
and hence, indicate a liquid vortex-like state. 

We first focus on the angular range (9 < 9c) where a di- 
vergent conductivity with vanishing current is observed, 
and we consider the scenario where the BG phase is sta- 
ble. Thus, one expects that excitations of some localized 
vortex lines lead to a nonlinear resistivity given by 1^ : 



p{J) = po exp{ - EK{Jc/J)''/kBT) 



(1) 



where po is a characteristic flux-flow resistivity and 
Ek{JcI J)^ represents the barriers against vortex motion. 
This expression is predicted to hold for various regimes 
of different behavior as the current density probes differ- 
ent length scales in the BG. Here, Ek acts as a scaling 
parameter for the pinning energy and Jc is the charac- 
teristic current scale of the creep process. When the cur- 
rent density is large enough that the growth of vortex- 
loops excitations of a line from its pinning track does 
not reach out the neighboring tracks, the half-loops ex- 
citations are relevant and lead to an exponent p = 1 
and Jc = Ji = U/{^od), where U is the mean pinning 
potential and d = y/^o/R^. With decreasing current 
density, critical size of half-loops increases with the re- 
sult that some disorder in the pinning potential becomes 
relevant. This situation yields the VRH process charac- 
terized by /i = 1/3 and by another important current 
scale Jc = Jo — 1/ {^og{p)d^), where g{p) denotes the 
density of pinning energies at the chemical potential of 
the vortex system. One expects that both Ji and Jq are 
insensitive to H±. We therefore consider that the only 
important effect of H± is to lower the scaling parame- 
ter for the pinning energy in Eq. (1), according to the 
formula: 



(2) 



where Ek is the mean kink energy, and the energy gain 
due to the tilt is obtained from the isotropic result 
^odH± by applying the scaling rule Eq. (3.12) of the 
Ref. §. 

We now present the following two-step method of fit- 
ting Eqs. (1) and (2) to the data, which allows to get 
a good understanding of the physics in our experiment. 
First, we plot in Fig. 2 the natural logarithm of V/I 
versus tanfl for / fixed. Note that tan9 happens to be 
here directly proportional to H±^ since / is maintained 
constant. A linear variation ln(y//) = A + Bta.n9 (ex- 
pression designated below as El) with A and B being 
current-dependent parameters is found for 9 varying up 
to 9c (solid lines), whereas above 9c data deviate from 



this behavior (dotted lines). This finding is another ar- 
gument supporting a true angular transition at 9c, as 
suggested in Fig. 1. Here, it should be noted that we 
draw such a conclusion from the observation of two in- 
dependent regimes of different behavior. One leads to a 
vanishing linear resistivity as 9 ^ 9^^ , and the other is 
evidenced by probing the regions of nonlinear resistivity 
above and below 9c- Another result is that B ~ with 
p undergoing a jump from 1/3 to 1 at /2 ~ 60 mA. Such 
a current crossover is clearly visible in the insert of figure 

1 /3 

2 where we plot together B normalized by /q and B 
normalized by Ii versus I depending upon whether / is 
< I2 01 > I2, respectively. Here, Jq and Ji are two cur- 
rents evaluated on the basis of the BG model, as will be 
seen below. It may be seen that below I2 , the solid line of 
slope ^ —1/3 fits the data very well, whereas the dashed 
line of slope ~ — 1 fits them best above l2- It therefrom 
follows that Bi/3{I) = k(/o/J)1/3 and Bi{I) = nli/I 
where the labels 1/3 and 1 differentiate between two fits 
below and above I2, respectively. In both these equa- 
tions, K w 0.6 is a dimensionless constant. It should 
be noted that such a crossover at I2 is clearly observed 
from V/I — I curves as well (see Fig. 1), where a sudden 
increase of V/I suggesting an increasing vortex motion, 
occurs at 12- Second, we plot in Fig. 3 and its insert the 
natural logarithm of V/I versus /"^/^ for I < I2 and 
versus I~^ for I > I2, respectively. We verify that the 
expression IniV/I) — C ~ D^I"'^ (expression designated 
below as E2) fits very well our experimental data (solid 
lines) with an exponent p ~ 1/3 for I < I2 and with 
^ = 1 for / > /2. In i?2, Dfj, is labeled using the above 
convention. The result is that C is a constant (within 
the experimental errors), while depends on 9. In Fig. 
4 and its lower insert, we show -D1/3 and Di versus tan^, 
respectively. A linear variation in tan^ is observed con- 
sistently with expression El. In upper insert of Fig. 4, 

1 /3 

we plot together Di/^/I^ and Di/Ii versus tan^. The 
result is that the data are superimposable onto a single 
straight line with slope « —0.6 consistent with the value 
of —K, as may be verified by identifying the expression 
El with the expression E2. 

In conclusion, V/I = Rq exp 
is clearly observed in our 
I2 (solid fines in Fig. 
Ro exp [- (k' - Ktan9) (h/I)] 
I2 (dashed lines in Fig. 1). In both fitting expressions 
K « 0.6 and k' « 0.4 are two parameters that a priori 
depend on the temperature and filling factor, Iq — 15 
A and Ii ~ 0.4 A are evaluated from BG model (see 
below), and Rq = e*^ « 1.1 pfl which is five orders of 
magnitude lower than the normal resistance. Therefore, 
for 9 < 9c we have rather strong evidence of two separate 
(albeit related) vortex creep processes peculiar to a BG: 
the half-loop expansion with p = 1 which is cut off by 
the crossover at I2 into the VRH process with p — 1/3. 
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is valid for data above 
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Moreover, we note that for each fiUing factor investigated 
in our experiment (/ = 0.13, 1/3, 2/3), we observe a 
continuation of the VRH process evidenced at 6* = 0° 
[ p2t . We therefore argue that the BG phase remains 
stable up to a critical tilting angle Oc, as predicted by 
NV |. 

To test the accuracy of the above view, we esti- 
mate from the theory 0,^, first the characteristic cur- 
rent scales, secondly the current crossover, and then 
the energy barriers in Eq. (1). Jq only depends on 
g{j2)- Although a form of g{x) is not yet available, 
an estimation of g{jl) can be done in terms of 7, the 
bandwidth of pinning energies due to the disorder, 
and hence, g{fL) « 1/ (^^7)- We take 7 = td + ji 
where ji arise from some on-site disorder, and t^ w 
U/y/Ek/kBTcxp {-V2Ek/kBT) is due to structural 
disorder ||^. Assuming the variation in the defect diam- 
eters as the first cause for random on-site energies, we 
approximate 7; to the width of the distribution of pinning 
energies P{Uk) = P{cK)dcK / dUx where the pinning en- 
ergy and the defect radius Ck are related to one another 

through the formula Uk = eoln^l + (c^/v^Cafc) 
and P{ck) designates a probability density function 
for the defect radii. In the following estimation of 7; 
, we use realistic Gaussian law centered at cq = 45 
A with a standard deviation of 6 A, as shown in fig- 
ure 1 of Ref. H. Then, taking Ek = d^/IJl with 
ei ~ eeo In (ao/^afc) where e is the anisotropic parame- 
ter and €0 = (1)1/ (AnnaXl^), and with U = Uof{T/T*) 
where T* — max [cq , V2£,ab] V^iUo is the energy scale 
for the pinning and f{x) = x^/2 exp (— 2a;^) accounts 
for thermal effects, and using appropriate parame- 
ters for BSCCO {Xab ~ 1850 A, ^ab « 20 Aand 
e « 1/200) we obtain Jq « 7/(^0^) « lOM/m^ and 
Ji ss U/{(j)od) « 2.6.10^A/to^. Assuming uniform cur- 
rents into the sample we have Jq ~ 15 A and /i « 0.4 
A. The basis for the estimate of the current crossover 
J2 between the VRH process and the half-loops regime 
corresponds to the strongly dispersive situation with 
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h/Io 



2.6.10- 



^ 1, so that 



Uh = Ji/Jo 

h ~ {U/iY'^h ~ {Uhf'^lQ with h < h < lo §, 
and thus, we obtain I2 ~ 6.10~^A in excellent agreement 
with the experiment (see Fig. 1). Below and above this 
current crossover, the energy barriers against vortex mo- 
tion are Ekilo/lY^^ and Ek{Ii/I), respectively, where 
Ek is given by Eq. (2). Using again the above usual 
parameters of BSCCO, we estimate Ek/ksT « 0.5 and 
e^cjiodHx/kBT « 0.7 tan 0, which correspond to the val- 
ues of our fitting parameters k' and k (see upper insert in 
Fig. 4). Thus, the experiment is also in good agreement 
with the energy barriers theoretically predicted below 
and above J2. 

Finally, we pass on to the linear regime observed 
above 9c (see Fig. 1). For 9 > 9c, Hwa et al. pre- 
dict on the basis of the CIC transition, the appear- 



ance of free moving chains of kinks oriented in the 
H± direction leading to a critical behavior of the lin- 
ear resistivity p ~ (tan 6* — tan 6'c)'^ characterized by an 
exponent = 1/2 (or = 3/2) in (l-fl) [or (2-fl)] 
dimensions. This type of behavior with v ~ 1/2 is 
evident in Fig. 5 which shows a plot of R/Rc versus 
tan6' — tan^c for three temperatures less than Tbg(O). 
Here, Rc is a scaling parameter comparable to Rq, and 
the critical tilt angle 9c{T) has been determined fitting 
R = R(j{ta,n9 - tan 6'c)^/^ to the data. In Fig. 2, the 
arrow indicates 9c obtained from this fit. The insert of 
the figure 5 shows Tbg{'^) — Tbg{0) ITbg{'^) as a function 
of tanfl where Tbg{(^) is obtained through the inversion 
of 9c and T. The solid line shown in insert is a fit to 
the data following Tbg(O) - Tbg{6) - \ i'An9\^/''^ with 
v±^ = 1.0 ± 0.1, as recently suggested by Lidmar and 
Wallin from numerical simulations |14| and observed in 
{K, Ba)BiO^ Note that such a value of vx^ is excel- 
lently consistent with the result previously found using 
the scaling theory of the BG transition at 6* = 0° jlj] . 

In summary, we have shown via a detailed quantitative 
analysis of the transport properties in BSSCO single crys- 
tal with columnar defects that a TME is supported by 
the experimental evidence of the creep processes (VRH 
and half-loop) showing the stability of the BG phase 
when H is tilted away from the column direction. Above 
the critical angle, the commensurate-incommensurate 
transition have been evidenced. Finally, our results in- 
dicate that the vortices behave as well-connected lines 
in accordance with the idea that columnar defects ef- 
fectively increase interlayer coupling in highly layered 
superconductors [ p^ , p7[ . 
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FIG. 1. V/I — I curves for tilted magnetic fields in an ir- 
radiated Bi2Sr2CaCu20s+y crystal. ^From the right to the 
left = 0, 5, 10, 15, 20, 25, 30, 35 and 40°. The curved 
lines are a fit of the Bose glass theory to the data. 
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FIG. 2. Angular dependence of the natural logaritmic 
of V/I for different currents. The solid lines are a fit of 

\n{V/I) = A-\-Bta.n6 to the data; the dotted lines arc a guide 
for the eye. Note that as the critical angle is approached from 
above, the linear resistivity drops continuously to zero (see 
Fig. 5). Insert: B/ll'^ vs. I (left axis) and B/h vs. I (right 
axis) with Iq and 7i determined from theory. The straight 
lines arc least-square fits whcrcfrom wc obtain /i = 0.31 ±0.05 
(solid line) for I < I2 and jj, = 1.0 ± 0.1 (dashed line) for 

I>l2. 
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FIG. 3. Natural logaritmic of V/I vs. I'^^^ for I< h and 
< 9c- Insert: natural logaritmic of V/I vs. I~^ for I > I2 

and 6 < 9c- In both plots the lines are a fit to the data using 
the form ln(V^//) = C - O^/"" with = 1/3 or 1 according 
to whether 7 is < /2 or > 72. 




Tane 

FIG. 4. Angular dependence of the coefficient defined 

in the regime I < I2 (see Fig. 3). Lower insert: angular de- 
pendence of the coefficient Di defined in the regime I > I2- 
Upper insert: angular dependence of Di/^/I^^^ and Di/Ii . 
In each plot the straight line is a least-square fit. 
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FIG. 5. Critical behavior of the linear resistance at dif- 
ferent temperatures for 9 > 9c{T). The curved line is a 
fit to the data using 7? = 7?c (tan^ — tan^c)^''^ in accor- 
dance with a CIC transition scenario. Insert: Angular de- 
pendence of the Bose-glass transition. The solid line is a fit 
of rBG(O) - Tbg{9) ~ I tan 61 1 1/"^^ with i^x = 1-0 ± 0.1 to the 
data. 



